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SIGNIFICANCE  AND  EXPLANATION 


Parallel  systems  arise  naturally  in  practice  in  engineering 
and  physics.  Therefore  it  is  of  substantial  significance  to  be 
able  to  efficiently  utilize  data  obtained  on  individual  components 
for  the  purpose  of  obtaining  an  overall  assessment  of  the  reli¬ 
ability  of  the  system.  The  methods  of  this  paper  may  be  employed 
for  this  purpose. 


Tii.  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


OPTIMAL  UPPER  CONFIDENCE  LIMITS  F''R  PRODUCTS  01’  P>  o  <;£./;  PYRAME"  ! 
WITH  APPLICATIONS  TO  Ti!H  INTERVAL  ESTIMATION 
OF  THE  FAILURE  PROBABILITY  OF  PARALLEL  SYSTEMS 


Bernard  Harris  and  Andr'-w  P.  Somr 


1.  Introduction  and  Summary 

A  problem  of  fundamental  interest  to  practitioners  in  relia¬ 
bility  is  the  statistical  estimation  of  the  reliability  of  a 
system  using  experimental  data  collected  on  subsystems.  In  this 
paper,  the  subsystem  data  available  consists  of  a  sequence  of 
Bernoulli  trials  in  which  a  "one"  is  recorded  if  the  subsystem 
functions  and  a  zero  is  recorded  if  the  subsystem  fails.  ?h 
for  each  of  the  k  subsystems  composing  the  system,  the  data  pro¬ 
vided  consists  of  the  pair  (n^,Y^) ,  i=l,2,...,k,  where  is 

binomially  distributed  (n.,p,).  We  assume  that  Y.  ,  Y  i .  a:-, 

mutually  independent  random  variables. 

The  magnitude  of  interest  in  this  problem  is  easily  evidv..o. 
by  the  extensive  literature  devoted  to  it.  In  this  regard,  se. 
the  survey  paper  by  Harris  (1977)  and  Section  10.4  of  the  bo ox  by 
Mann,  Schafer,  and  Singpurwalla  (1974).  In  addition,  the  Lofeuso 
Advanced  Research  Projects  Agency  has  recently  issued  a  Handbc 
for  the  Calculation  of  Lower  Statistical  Confidence  Bounds  on 
System  Reliability  (1980) . 

Historically,  the  first  significant  work  on  this  problem  wu 
produced  by  Buehler  (1957).  However,  Buehler’s  method  os  des¬ 
cribed  in  that  paper  is  difficult  to  implement  compulaticra ■ 
when  k>2. 
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In  this  paper,  we  examine  the  problem  of  obtaining  upper 


confidence  limits  for  products  of  Poisson  parameters.  This 
problem  is  studied  by  means  of  majorization  methods  and  Schur- 
convexity,  such  as  described  in  the  book  by  Marshall  and  Oikin 
(1979) .  A  significant  application  is  the  determination  of  con 

dence  limits  for  the  reliability  of  systems  of  k  parallel  sub¬ 
systems,  a  fundamental  problem  in  the  statistical  analysis  of 
reliability. 


i 


3 


2.  Exact  Solutions  for  Products  of 
Poisson  Parameters  for  Small  Failure  Combinations 

Let  X  =  (X^,X^ , . . . ,X^)  be  independent  Poisson  random  variables 

k 

with  parameters  A,  ,  ,  . . .  ,  X,  ,  k  >  2,  and  let  h(A)  =  IT  A..  Let 

1  ^  K  -  i=l  1 

k 

g(x)  =  H  (x.+d)  ,  l<d<1.5  ,  x.=0,l,...  (2.1) 

i=l  1  1 

and  denote  the  ordered  points  in  the  range  of  g(x)  by 
j i <3 2 < • • • <3m< • • • •  Define 

Aj_  =  |x  | g  (x)  =  j±|  .  (2.2) 

Since  x^,  i=l , 2 , . . . , k,  takes  on  non-negative  integral  values, 
v?  regard  it  as  desirable  to  have  d  in  (2.1)  only  assume  non- 
.ir.teger  values.  This  has  the  effect  of  making  the  partition 
defined  in  (2.2)  finer  than  would  be  the  case  if  d  were  an  integer. 
It  is  easily  verified  that 


q  =  sup-,  .  q .  (  i 

*■  i=l  1  xi  Aj_ '  ■‘-in  j 


x  a  .  n  -x  x .  . 
j  i  ]  i  13 


X  -1 


J  .  J  =  a 
- 1 


Thus  in  practice  a  /  .'!  n.  may  be  employed  as  an  approximate  (1-r 

n  i=l  1  k 

upper  confidence  limit  for  q.,  q.=l-p..  In  this  sense  the 

i=l  1  1  1 

methods  of  this  paper  can  be  used  as  approximations  for  estjmati 
the  reliaoiiity  of  parallel  systems  when  independent  bmomially 
distributed  data  is  obtained  for  each  component.. 

We  proceed  by  showing  that  a • x)  is  a  Schur-concave  function 


and  consequently 


|x  |  g  {>:)  _  g  (  xq  )  j 


is  a  Schur-convex  set  (see  Marshall  and  01  kin  (1974),  pp.  1139-9 
and  Nevius,  Proschan  and  Sethuraman  (1977),  p.  264).  The  Schur- 
concavity  of  g(x)  follows  immediately  by  noting  that 


..  .  f  9g(x)  9g(x) 1  , 

Xn  J  - -  — - I  " 


kl  2  !  3x 


5x  j  - 


<  0  . 


Define  F(x  ;A)  by 
o 


F  ( x  ;  \  )  ~  T  f  ( x  .  ;  A )  =  P-(b-  ) 

°  x . G  B  ~  1  A  XC 

1  xo 


and  let 


u(x  ;a)  -  sup  F  (x  ;>.)  ,  0<a<l  .  (2. 

0  h ( X ) -a  ° 

Since  the  Poisson  distribution  has  a  monotone  likelihood  ratio, 

u(x  ;a)  is  a  strictly  decreasing  function  of  a  for  fixed  x  . 
o  J  o 

Hence  for  every  c,  0<c-'l,  there  is  a  unique  a(c)  such  that 

u  ( x  ;  a  ( c )  )  =  c  .  ( 2  .  ' 


Consequently,  we  also  have  that  a  (sec  (2.3))  is  the  solution  i 


a  of 


u(xQ;a)  -  a  . 


'2. 


(2.8)  is  established  exactly  as  in  Harris  and  Soms  (1330'. 

The  method  ... ogy  to  be  employed  is  as  follows.  If  Fix;  •.!  is 

a  Schur-concave  function  of  =  -In  A^ ,  i=l,2,...,k,  then  it 

1/k- 

follows  that  u(xQ;a)  =  F(xo;a  1)  ,  where  1  -  (1,1,...,1),  a no 

then  the  solution  in  a  of  u(x  ;a)  =  a  is  an  optimal  upper  cor. fi¬ 
le 

dence  limit  for  n  A..  This  will  entail  verifying  (for  fixed  x  } 
i«l  1 

that 


<VV 


liF(y.:Xt  3F  ( X  ;  ' )  ) 
|  _ 0_ _ _  _ Q  | 

R,  *  3  R. 


<  0 


(2.9: 


^  °  "1  °  iX2  J  - 

(see  Marshall  and  Olkin  (19  74)  ,  p.  1190).  Accordingly  wo  hat 
following  theorem. 


k 

Theorem  2.1:  Let  g(x)  =  II  (x.+d),  l<d<1.5,  k>3.  Define  0_,  us 

i=l  1 

the  j-vsetor  all  of  whose  components  are  zeros.  Then  let 

x(1)  =  0k,  x(2)  =  (l,0k_1),  x(3)  =  (2,0k_1),  x(4)  =  (l,I,0k_2  • 

x(3)  =  (3,0,,_1),  x(6)  =  (4,0k_1),  x(7)  =  (2,l,0k_2)  and 

x<8)  =  (5,0.,  .).  The  set  A.  defined  by  (2.2)  is  the  point  x ' 1  ’’ 

and  the  different  permutations  of  its  components,  i=l , 2 . 

Further,  for  j=l,2,...,7,  F(x*'^;A)  is  Schur-concave  in  , 


i-1, 2, 


,k. 


Proof :  In  the  sense  of  the  ordering  given  by  (2.2),  obviously 

x(1)  <  x(2)  <  x(3)  <  x(5)  <  x(6)  <  x(8).  Trivially,  d  (2+d)  <  Q-td) 

and  hence  (2+d)dK_1  =  g(x(3))  <  (l+d)2dk_2  =  g(x(4)  )  .  Similar  ly  , 

2  -  (4)  ~  (S' 

since  l<d<1.5,  (1+d)  <  d(3+d)  and  hence  g(xv  ;)  <  g(xv  ') 


In 


the  same  way  g(x^)  <  g(x^),  g(x^)  <  g(x^),  g(x^)  <  g(  1,1, 1,0. 


and  a  ( x 


(8) 


g(2,2,0k  2)  '  establishing  the  first  part  cf  the  cone  1 


i 


fei 


6 


In  order  to  establish  Schur-concavity ,  v/c*  must  verify  (2.9) 
Thus  consider 


k  -X  .  xf^/  ... 

I  n  e  1  A.  '  /  x.™ 

.£A.,j<«.  i=l  X  /  1 


3  3 


[  2  .  i  j ) 


“R  * 

where  =  e  1.  Define 


k  -R. 


G(x(£);R)  = 


a  F  ( X  '  ^  '  ;  a  )  _  3F(X(£)  ;>.))  /  i=l 


-  y 


3R, 


3Ro  J 


r  2 . 1 1 ) 


Letting  R  =  (Rx , . . . , R  ) ,  we  obtain 


-  (1)  ~  '^1  -R-> 
G(xu;;R)  =  (e  1  -  e  2)  , 


)  9  \  —  — Ri  ~ R k  -R . 

G(xu;  ;R)  =  (e  1  -  e  2)  (  l  e  x)  , 

i=l 


/  ^  \  * R -R_  k  -R .  k  2Ri 

G  (x 1 J  ;R)  =  (e  1  -  e  2)  (  y  e  1  +  J  e 


i=3 


i=l 


)  , 


~  (  4 )  -  *"?  ^  _  j.  is- 

G  (x  4  ;  R)  =  (e  1  -  e  2)  (  [  ®__  +  [  e  1  3 )  , 


»OD 

~ ^ k  _  i  k  -R .  -R . 


i=l 


K] 


_  p  — R  v  “2R.  R_—R_ 

_  { 5 )  Ki  0  l  ^  1  2 

G(xlb;  ;R)  =  (e  1  '  e  2)  {  l  +  ®__ - 

i=3  ^  2 


k  -R . -R  k  -2^i 


+  I  e 

i< j*  (i, j)*(l,2) 


i  3  +  r  e 
i=l 


+  ^  3!  }  ' 


_r  j.  _2Ri  -R. -R_ 

G(ix6)  ;R)  =  (e  1  -  e  2)  (  l  -  +  § - - - 

i=3  2  2 


X  -R.-R. 

v  il 

L  e 

i<  j  »  (i,  j)?«(l,2) 


k  _3Ri  -2RrR2  -Rr2R2 


4  ) 

i  =  3 


+  e 


3! 


3  ! 


.  -4R. 

*  e  1 

+  1  ) 


4  ! 


and 


7 


.  -R:  -R2  n"2Rl'R2  /R1'2R2 

G  ( x  '  '  ;R)  *  (e  -  e  l)  (- - ^ -  +  ~3 - + 


-R .  -R 
,  1  3 


3<  l-'  j 


k 

V 


-2R . -R . 
3 


(i,j)^(l,2)or(2,l) 


-R  -R  -R,  -R2 

N’ow  R., '  R0  implies  e  <  e  and  thus  (R^-R2)  and  (e  -e  ) 

have  opposite  signs.  Hence  it  follows  that  FCx^1^;!),  i=l,2,3, 

4,5,"  i5  Schur-concave  in  R..  The  verification  that  F(x^;A)  is 

1  -R2  ~Ri 

Scnur-cor.oave  nay  be  accomplished  by  letting  k=2,  e  =  ce  , 

c>0,  and  examining  the  discriminant. 

"i~  shew  that  (2.9)  need  not  be  positive  for  all  xq,  consider 

k  ■-  2  and  x  =  (7,0).  Then 


'  R  i  Rp  e 

G ( x  ;R)  =  (o  1  -  e  A)  (- 


-6R,  -6R, 

_ +  e 

6! 


-4R. -R-  ~3R. -2R_  — 2R^-3R2  — Rj  ^^2 

^  4-  ^ _ +  e _  _ 

~  “  31 

-3R.-R,  -2R  -2R?  -R,-3R2  -2R1-R2  -r1-2R2 

2 _ _ _ lt£ _ i® _ +  s - p - ) 

4!  3 

“R1  _R2  . 

and  th  is  is  Scrhur-ccnvex  near  e  =  e  “4. 

3  .n-hler  provided  an  extensive  discussion  of  this  problem  for 

th;  ordering  function  determined  by  the  product  of  the  upper 


confluence  limit: 


for  the  individual  components.  In  particular, 


:e  ~v  ided  some  numerical  tabulations  for  k=2.  Asymptotically 
hr  •hicr ’ s  ordering  function  is  given  by 


g  (x)  =  H  (x.+z  ,x/  )  , 

B  ,=1  i  a  i 


i-a-u7  ,  z  satisfies  $  (zq)  =  1-a  ana  l(x)  is 
;  normal,  cumulative  distribution  function.  It  is  easy 
is  Schur-concave  (see,  e.g.,  Marshall  and  Olkin 


3.  Bounds  on  Confidence  Limits 


In  this  section  we  employ  major ization  techniques  decor  jo 
in  Proschan  and  Sethuraman  (1977)  and  Nevius,  Proschan  and 
Sothuruman  (1977)  to  obtain  bounds  for  a  .  Throughout  this  s-/ 
••.Los  wo  assume  only  that  the  ordering  function  g(x)  is  strio.  1 
increasing  in  each  component  and  Schur-concave  and  thus  the  r.o 

B -  will  be  Schur-convex  (see  the  discussion  immediately  preco 

o 

inu  (2.5))  . 

In  order  to  proceed,  we  need  the  preliminary  results  esta 
lished  below. 

1/k 

Theorem  3.1:  Let  c  and  a  be  given  with  oka  '  and  consider  t 
set  A(a,c)  of  vectors  A  =  , A^) ,  A^>0,  such  that 


n  A.  =  a 
i=l  1 


and  j-  A .  =  c 
i-1  1 


fL  * 

Let  S.  =  _  max  l  A..  Then  there  is  a  unique  A  sA(a,c)  of 
^  AeA(a,c)  i=l  1 

the  form  A.  =  M.,  l<i<j,  A.  =  m.,  j+l<i<k,  M .  >m .  ,  S.  =  j.M .  . 

i  3  -  ~J  i  3  “  "  3  3  3  J  3 

1/k 

Proof :  The  condition  oka  '  is  a  consequence  of  the  arithmet 

geometric  mean  inequality  and  insures  that  A(a,c)  is  non-triv? 
for  k>3.  If  k=2,  there  is  only  one  solution  of  (3.1)  with  A N 
and  hence  the  Theorem  is  trivially  true.  Consequently,  suppos 
k>3.  Then  for  fixed  j,  (3.1)  requires  that  any  solution  of  th 
required  type  satisfy 


[ 

I 


f  j  (M)  =  [  (c- jM)  /  Ik-  j  )  ]  k  \  l^j<k-l,  O^M-'c/  j  . 


Note  that  f.  .  (0)  =  f.(c/j)  =  0,  and 
j  3 


(3.2) 


...  k-j-1 


f.(M)  =  (c-MK)(^— -)  (~±- -) 


Thus,  f j (M)  is  increasing  for  0<M<c/k  and  decreasing  otherwise, 

further  f^  (c/k)  =  (c/k)  ^>a.  Hence  there  is  ^exactly  one  solution 

M.  of  f.(M)=a  with  M.>c/k,  and  therefore  M . >m . . 

3  3  3  3  3 


(3.3) 


Now  assume  that  for  some  j,  l<j<k-l,  the  vector 

*  *  *  *  3  * 

A  =  ( A.  ,  A_  ,  .  .  .  ,  A,  )  with  I  A.  =  S.  is  not  of  the  foi 
lx  K  i  3 


s  not  of  the  form 


(M . , . . . ,M . ,  m . , . . . ,m . ) .  Then  let  A, .  =  S ./ j  and  A_ .  =  (c-S  ) / (k- ' ) 
3  3  3  3  A3  3  <-3  3 

" - v - "  v- - ^ - / 

j  k-j 

Define  i!  -  . Ajk>  by  it  -  q.  ,  l<l<j,  i’j  -  ~2j  , 

j+l<i<k.  Since  the  geometric  mean  of  a  set  of  positive  numbers 

whose  sum  is  fixed  is  a  maximum  when  they  are  all  equal,  we  have 

k  ,  . 

HA..  > a.  Now  A.  is  of  the  required  form,  however,  from  (3.2)  and 
i=l a  1  k  ,  3 

(3.3),  II  A..  >a  implies  that  there  is  another  solution  of  the 
i=l  1,1 

required  form  with  A ^  >  S^/j,  l<i<j,  contradicting  the  maximal ity 

of  S  .  . 

3 

From  (2.5)  and  (2.6),  we  can  write 


u(x  ;a)  =  sup  P~ (B~  )  =  sup 

O  1  A  X  , 

k  o  c  k 


Pr  B~  .  3.4) 

A  X 

O 


n  A  .  =  a 
i=l  X 


y  A . =c , n A . =a 

■  -i  i  i 

1=1  1=1 


We  state  now  the  main  result  of  this  section,  using  Theorem  3.1. 

Theorem  3.2:  Let  v.  =  M. ,  v.  =  iM.-(i-l)M.  2<i<k-l, 

~~  k-1  1  11  1  1-1 

v.  =  c  -  \  v.,  where  M.  is  specified  by  Theorem  3.1.  Then 

K  i=i  1  1  - 


]  0 


u (x  ;a)  <  sup  P~ (B-  )  . 

C  V  X 


(3.5) 


a 

Proof :  Since  ) 

i=l 


k 

l 

i=l 


v.  =  S  .  ,  l<j<k-l,  l  v.  =  c,  v  majorizes  every 
1  J  —  “  -•  -  ■>  1- 


k  k 

A  with  1  >  ■  =  c ,  H  X.  =  a  (Theorem  3.1).  Then  (3.5)  follows, 

i=l  1  i=l  1 

since  if  A^  majorizes  ,  then  for  any  Schur-convex  set  A, 

Pr  (A)  >  P“  (A)  (Proschan  and  Sethuraman  (1977)  and  Nevius, 

A1  ~  a2 

Proschan  and  Sethuraman  (1971),  p.  264  and  pp.  267-9). 

The  vector  v  may  be  interpreted  as  the  best  vector  that  major 

k  k 

izes  all  vectors  A  such  that  7  A.  =  c  and  H  A.  =  a.  More 

i=l  ^  __  i=l 

specifically,  there  is  no  vector  w  ^  v  such  that  v  majorizes  w  and 
w  majorizes  all  A  satisfying  the  two  conditions  given  above. 

The  following  is  a  suggested  method  for  employing  Theorem  3.2 
Find  a.  such  that 


a 


F  (x. 


1A 


1) 


Next  calculate  the  smallest  a,  say  a^,  such  that  sup  P~  (B^  )  <  ct . 

c  v  xo 

If  a  =  a,,  this  is  the  exact  solution.  Otherwise  a ,  < a  and 
m  d  dm 

suo  P  (B-  )  <  a  (here  a  =  a  )  and  the  solution  a  satisfies 
v  x  m  n 

c  o 

a(jlanlam'  The  vector  v  may  be  calculated  by  any  of  a  variety  of 
numerical  techniques.  In  the  numerical  examples  presented  here, 


interval  bisection  was  employed. 


Example  1 :  Let  k  =  5,  a  =  25,  c  =  15.  Then  the  4  vectors  A^,  A2/ 

A-  and  A.  of  Theorem  3.1  are 
3  4 


1 

"l 


(9.9660, 

1.2585, 

1.2585, 

1.2585,  1.2585) 

(6.2004, 

6.2004 , 

.  8664  , 

.8664,  .8664) 

(4.6696, 

4.6696, 

4.6696, 

.4955,  .4955) 

(3.7172, 

3.7172  , 

3.7)12, 

3.7172,  .1309)  , 

'c 


and 
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from  which  v  is  determined  to  be 

v  =  (9.9660,  2.4349,  1.6079,  .8601,  .1309)  . 

Note  that  in  the  above  example  Vj_>v^>. .  .^v^.  This  in  fact 
is  always  true,  as  the  following  theorem  establishes. 

Theorem  3.3:  For  v  defined  by  Theorem  3.2,  we  have  v^^v2- •  •  •  >vjc* 

Proof :  It  follows  immediately  that  — V2 '  s^nce  ^1-M2"  Consider 

therefore  v ^ ,  j>2.  Vj>v^+^,  j= 2, 3, . . . ,k-l  holds  if  and  only  if 

>  ( j+l)Mj+1~jMj 


jMj  >  ((j+l)Mj+1+(j-l)Mj_1)/2  , 

where  Mk  =  c/k  (satisfying  the  condition  =  c  =  kM^  of  Theorem 
3.1)  . 

Let  XA  =  (1— a j )  A  j=  2,  3, . . .  ,k-l,  where 

Ct  j  =  ( 1/2 )  +^(  1/  (2  j ) )  and 

Aj  =  (Xj1,Xj2,-.  • 


X..  =  M.  ,  l<i<j,  X..  =  m.  ,  j+l<i<) 

:i  ]  - J  3i  3  J  — 

It  follows  that 


n  X  i  >  a  , 
i=l  Aj' 

k 

since  \  In  x.  is  a  concave  function  of  x.,...,x..  Now  let 
i=l  •  1  k 

AB,,i  =  XA i=1'2 . 3'  XByi  =  (i4+1XAj,]/<k“j)  ’ 

i=j+l, . . . ,k.  Then 

3  j 

rc  x„  i  >  n  Xj.  i  ,  j=l , 2 ,  .  .  .  , k-l  . 

i=l  Bjrl  i=l  Ajfl 


Thus,  using  the  properties  of  in  Theorem  3.1, 

jM j  >  (j-1)  f  +  (l-aj)m^_1  +  a.jMj+1  , 

yielding 

jM ^  >  ( (j+l)Mj+1+(j-l)M._1)/2  +  (1-om )mj_1  , 
which  establishes  the  theorem. 

To  illustrate  the  techniques  of  this  paper,  we  compare  numer¬ 
ical  values  obtained  by  the  above  method  with  those  given  in  the 
examples  from  Mann,  Schafer  and  Singpurwalla  (1974,  p.  505).  Fro; 
now  on  we  assume  d  =  1.1. 

Example  2 ;  For  xq  =  (1,2,1)  we  obtain  a^  =  =  20.56  for  a  =  .10 

In  Mann,  Schafer  and  Singpurwalla,  an  AO  non-randomized  confidence 
bound  of  20.7  is  obtained. 

Example  3 :  Let  xq  =  (2,3,5),  a  -  .10.  Then  we  obtain  a^  =  135.46 
A  summary  of  computer  calculations  which  establishes  135.46  <  < 

142.46  is  given  below  in  Table  1.  With  the  exception  of  the 
likelihood-ratio  value  of  133  and  the  AO  non-randomized  confidence 
bound  of  129,  all  the  other  confidence  bounds  given  in  Mann, 
Schafer  and  Singpurwalla  exceed  the  upper  bound  of  142.46.  For 
k-3  it  is  possible  to  do  a  direct  computer  tabulation  of  u(x  ;a). 
This  gives  an  =  135.46,  the  diagonal  value. 

Insert  Table  1  here. 

The  two  examples  below  are  for  four  and  five  component  system 
for  which  there  are  no  comparable  numerical  examples  available. 

Example  4:  Let  x  =  (2, 2, 2, 2)  and  a  =  .10.  Then  a.  =  a  = 

°  an 


150.63 


13 


Example  5:  Let  xq  =  (2, 2, 2,2, 2)  and  u  =  .10.  Then  a^  =  429.69. 

A  summary  of  the  computer  calculations  which  establish  429.69  ^ 
an  <  435.69  is  provided  in  Table  2. 

Insert  Table  2  here. 

1/k 

As  a,  increases,  the  difference  between  a,  and  a  becomes  wide1' 
d  d  n 

Thus  the  techniques  of  Section  3  are  more  useful  for  small  x^,or  equiva- 

1/k 

lently,  small  a^  '  .  For  example,  for  xq  =  (5,5,5),  a^  =  387.18, 

and  it  is  not  practical  to  compute  a^  because  it  is  much  bigger 

than  a,.  However,  direct  tabulation  of  u(x  ;a)  reveals  once  more 
d  o 

that  a^  =  u(x  ;a).  A  justification  of  why  a^  =  u(xo;a)  for  large 

1/k 

a^  is  given  in  the  Appendix.  This,  together  with  the  results 
of  Section  2,  suggests  very  strongly  that  for  all  practical 
purposes  a^  =  an> 

Remarks:  Note  that  Tables  1  and  2  are  virtually  linear  in  their 

behavior  in  the  neighborhood  of  the  solution.  This  suggests  that 
solutions  are  obtainable  by  interpolation  and  then  one  should 
subject  them  to  verification. 

The  calculations  described  above  utilized  two  short  FORTRAN 
programs  for  2-10  components.  Listings  are  obtainable  from  the 
authors . 


4.  Comparisons  with  Buehler's  Tables 

In  order  to  provide  an  illustration  of  the  performance  of 
k 

g(xQ)  =  n  (x^+d) ,  l<d<1.5,  when  compared  with  the  tables  given 
by  Buehler  (1957),  we  chose  d=l.l,  k=2.  For  k=2,  the  values  of 
an  and  a^  coincided  for  both  the  ordering  based  on  g (x)  and 
Buehler's  ordering  and  further  were  for  all  practical  purposes 


equal  for  the  two  different  orderings. 


In  Table  3  wc  give  Buehler's  upper  confidence  limit,  Buehler's 

diagonal  value  and  the  exact  upper  confidence  limit  and  diagonal 

value  corresponding  to  g,  denoting  them  by  a  r  ,  a  , r  ,  a  and  a,  , 

J  nB  dB'  ng  dg 

respectively.  These  values  are  provided  for  all  failure  combina¬ 
tions  from  (0,0)  to  (5,5)  for  n=.l. 


Insert  Table  3  here 


An  examination  of  Table  3  shows  that  differences  between  the 
four  alternatives  presented  are  small  for  the  specific  example 
(k=2,  a-. 1) . 


5.  Concluding  Remarks 

In  this  paper  a  procedure  for  obtaining  bounds  on  an  optimal 
upper  confidence  limit  for  the  failure  probability  of  a  parallel 
system  is  given.  The  procedure  employs  the  theory  of  major- 
ization  and  xs  valid  for  an  arbitrary  number  of  components  and 
gives  the  exact  answer  or  narrow  bounds  when  the  observed  number 
of  failures  is  small  for  each  component.  In  addition,  numerical 
and  asymptotic  justif ication  is  given  for  using  as  an 
approximation  to  an«  Tables  of  a^  are  in  preparation  for  moderate 
numbers  of  failures  for  3,  4  and  5  components  and  will  be  avail¬ 
able  in  the  near  future. 
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Appendix 


Theorem  Al;  Let  X^,  1  _<  i  £  k,  be  independent  identically 

distributed  normal  random  variables  with  means  X  and  variances  X 

Let  X >  1  <_  i  <_  k,  be  independent  normally  distributed  random 

variables  with  means  t.  and  variances  t.,  where  t.  =  X  +  0(XC), 

i  i  i 


k 

c  <  l,as  X  -*■  00 ,  and  II  x  .  =  X  .  Let  &  be  given, 

i=l  1 

let  a  be  a  specified  positive  real  number,  let  Z ^ 
k 

Z_  =  n  (X0 .  +  a)  and  let  d(X)  satisfy 
*  j=l 

P[Z1  <  d  (X ) ]  =  6. 


0  <  &  <  1, 
k 

■  n  (X  .  +  a)  , 

j-1  3 


(A.  1) 


Then  as  X  -*■  °°, 


B  -  P[Z2  <  d (X)  ] 


0[  (£.n  X)1,5X-1]  ,  c  <.  0, 


Otx^1] 


(A. 2) 


,0  <  c  <  1, 


Proof ;  Throughout,  let  $  and  $  denote  the  density  and 
distribution  function  of  the  standard  normal.  Clearly, 


P[ZX  <  d (X )  ]  -  P[Z2 


where  x  =  (x^,x2#. 
of  xn'xi2'  *  *  *  'Xlk 


00  CO  00  ^  ^  ^ 

<  d(X>]  =  /  J  •••  /  (f.(x)  -  f 2 ( x ) ) dx , 

—  00  —00  —00 


•  'xk} 
and 


(A. 3) 


k 

(x:  II  (x .  +  a)  <  d  (X)  } 

j-1  3 

,  f^  is  the  probability  density  function 
f2  is  that  of  X21,X22 , . . . ,X2k.  Now 


16 


PtXlj  1  -a,  3 

and 

HX2.  >  -a,  j 

Consequently,  for 
m  >  0  such  that 

P[X. 

: 

Then,  for  i  =  1,2, 
P[Zi  <  d (X) ]  = 

and  therefore 

P[Z.  <  d  (X )  ]  -  P[Z. 

Next,  we  calculate 
?[Z1  <  d (X) ,Xlj  >  - 
Now 

p[Zl  <  d(x),Xlj  >  - 


-  l 


■2 . k!  -  i1  '  iWr  1,<^I72)) 

■f1/2 

1  •  r^r  ♦( rva 


—  1,2,<«  *,k]  >  IT 

"  j=l 


T  . 

3 


(A. 4) 


(A. 5) 


X  sufficiently  large,  there  exists  a  constant 


LD 


>  -a,  j  =  l,2,...,k]  >  1  -  e~mX  ,  i  =  1,2.  (A. 6) 


PtZ±  <  d (X) #Xij  >  -a,  j  =  1,2,.. ,,k] 


+  P[Z  <  d  (X ) ,  U  (X.  .  <  -a)] , 

J-  -.ill 


j-1  13 


<  d{X),X..  >  -a,  j  =  1 , 2 , . . . ,k]  <  e"mX.  (A. 7) 


■a, 


3  ~  1/2, • . . , k]  “  P  [  Z j  ^  d (X ) i%2j  ~  — a,  3  ~  1/2,..., 


■a,  j  1,2, . . . ,k| X^j  =  Xj ,  j  2 , 3 , . . . ,k] 


=  P 


X  <  d(X)  _ 

XU  -  ~  a 

II  (x.+a) 

j=2  D 


(A. 8) 


=  $ 


d(X) 


IT  (x.+a) 
j~2  ^ 


-  a  -  a 


/X 


1/2 


*(b~) 

X 


;V 
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Therefore 


P[Z1  <  d (X ) #X1 j  >  -a,  j  =  1,2, ...,k] 


=  /  /  $ {b~ ) g1 (x2 ,x3 , —  ,xk)dx2dx3  .  .  .dx^  , 

—a  —a  —a  1 


where  (x2 ,x3 , . . . ,x^)  is  the  probability  density  function  of 
X12 ,X13 ' ' ' * ,Xlk*  From  (A. 6),  we  have  that 


(A. 9) 


/  J  •••/  <f  (b~)g1(x2,x3, . .  .  ,xk)dx2dx 


_oo  _ oo  —00 


3 ' * *dXk 


(A. 10) 


“"ITlX 

-  /  /  •••/  <f  (b-Jg^^  (x2  ,x3  , .  .  .  ,xk)dx2dx3.  .  .dxk  _<  e“  . 


-a  -a  -a 


Hence  we  will  estimate  the  first  expression  on  the  left  hand  side  of 
(A. 10).  Similarly,  for  Z^  we  will  consider 


/  J  **•/  * 

— 0O  —00  —00 


d(A) 


II  (X.+T.) 


-  •  -TjAi72 


j*2 


3  3 


g2 (x2,x3, . . . ,xk)dx2dx3. . .dx^  , 

(A. 11) 


where  g2 (x2,x3» . . . ,xk)  is  the  probability  density  function  of 
X22 ,X23 '  *  •  *  ,X2k’  In  t*ie  ^^rst  integral  in  (A. 10),  let 
(y±  -  A)/A1//2  =  ui  and  in  (A. 11)  let  (yi  -  T /tV2  =  uif 

i  =  2,3,...,k,  obtaining 


18 


I  f 

—  00  —  oo 


-  /  I 


M  M 

/  ; 

-M  -M 


_  a  .  X  /xl/2 


\'  H  (x.+a) 

ir  3 


gl (x2 ,X3 ' * • *  »xk)dx2dx3 • ' *dxk 


•  i 


d  (X) 


k 

H  (x.+a) 


-  a  -  t. 


A 


1/2 


j  =  2 


g2 (x^ / * • • • ,x^) dx2dx^ . . . dx^ 


\  '* 
i 


d  (X ) 


H  (X1/2x  +X+a) 


-  a  -  X 


A 


1/2 


(A. 13) 


j  =  2 


d  (  \  ) 


-  a  -  t  . 


I  tt  r-1  /2x..+i^+a) 


A 


1/2 


3=2 


:>  j 


I  k 


n  4>  (x-i ) 

j  =  2  3  , 


dx„dx_.  .  .  .  dx, 
2  3  k 


+  V 

2 

1/2  fk-1)©-^  ^  -1 

where  M  =  (?.  In  X)  '  and  RM  <  4  . -•<  y_ -  =  0(X  A)  . 

M  "  <2,)1/2M 

Using  d(X)  =  xk  -  kd(X)Xk_1/2  ,  k  (X)  =  0(1), 


d  (X) 


-  -  a  -  X 


3 


1  /"> 

n  (X  x..+X+a) 

=  o  i 


/X1/2  =  (X1/2-k,  (X)  )  n  (l+x.X_1/2+aX_1)“1 
I  j  =  2  D 

-X1/2  -  ax'1/2  . 


Since  |x.  <  M,  we  have 

l  — 


(1 


+  x.x“1/2+aX-1)-1  =  l-x.X_1/2+(-a+x2)X“1+0(  (Jin  X)1,5X"1,5). 


.  A- 


Thus 


nl/2  .  ....  *  ^  ,“1/2  A  ,-1,-1  ,1/2  ,-1/2 

(X  -  k.(X))  II  (1  +  x.X  +  aX  )  -  X  -  aX 

a  j=2  11 

=  -  I  x.-k  (X)-kaX'1/2+k,(X) (  J  xi)X‘1/2+(^  x2)X~1/2 
1=2  1=2  1=2 

+  (  l  l  x . x  .  )X  ~1/2  +  0 ( (£n  X)1-5X'1)  . 
2<i<j  1  3 

Similarly,  using  x.  =  Xk/  II  x  .  ,  x./X  =  1  +  0(XC  1)  , 

1  j  =  2  3  3 

(x./X)1/2  =  1  +  0(XC_1),  j  =  1,2,. ..,k,  |xi|  <_  M,  we  have 


(A. 14) 


-k - ^ - a  -  x.  /x1/2 

i  n  (x1/2x.+x .+a) 

j=2  i  33 

=  ( x 2^/X x ) 1/2  (X  1/2-k^  (X ) )  (1+x  jX~1/2+axT1)  -1-X 1/2  -ax"1/2 


(A. 15) 


=  -  l  x . -k . (X ) -kaX-1/2+k . (X ) (  £  x.)\~x/*+(  [  x* 
i-2  1  d  Q  ^Z-2  1  i=2  1 


“1/2, 


2, ,-1/2 


+  (  ll  x.x.)X_1/2  +  0(XC"1)+0( (£n  X)1,5X_1). 
2<i<  j  1  3 


Combining  (A. 14)  and  (A. 15)  with  (A. 7),  (A. 9),  (A. 10)  and  (A. 11) 

establishes  the  theorem. 

For  c  <  J  standard  weak  convergence  arguments  show  that 


Um(8  -  P[Z?  <  d{X)] )  =  0. 

X-H» 


In  this  case  Theorem  A1  provides  additional  information  by  specifying 
the  rate  of  convergence. 
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By  standardizing  the  first  expression  in  (A. 13)  and  applying 
the  dominated  convergence  theorem  the  following  result  can  be 
obtained . 

Theorem  A2;  Let  x^i'  ^  ^  be  independent  identically 

distributed  normal  random  variables  with  means  X  and  variances  X. 
Let  X2i'  1  i.  i  £  k»  be  independent  normally  distributed  random 

c 

variables  with  means  and  variances  t^,  where  =  X  +  0(X  ), 

c  <  1  and  let  g,  and  be  specified  as  in  Theorem  Al. 

Then 

Um(6  -  P[Z,  <  d (X ))  =  0. 

X-+°° 
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1.  Summary  of 

Calculations 

Used  to 

Obtain  the  Upper 

Bound  for  a_ 

in  Example 

_3 

n  - 

a 

V 

sup 

P-  (B- 

t  r  V 

c 

V  A 

1 

135.46 

13.0680 

4.7283 

1.7108 

.1101 

136.46 

13.0867 

4 . 7409 

1.7173 

.1086 

137.46 

13.1053 

4 . 7532 

1.7240 

.1071 

138.46 

13.1239 

4.7656 

1.7305 

.1057 

139.46 

13.1423 

4.7780 

1.7370 

.1042 

140.46 

13.1607 

4.7902 

1.7435 

.1028 

141.46 

13.1789 

4.8024 

1.7500 

.1014 

142.46 

13.3299 

4.8057 

1.7325 

.  9999 
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Summary  of  Calculations  Used  to 
Obtain  the  Upper  Bound  for  in  Example  5 


a  sup  P~ (B~  ) 

-  ^  v  x 


429.69 

.1016 

430.69 

.1013 

431.69 

.1010 

432.69 

.1007 

433.69 

.1004 

434.69 

.1001 

435.69 

.0998 

1 
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-nB 

3.  Comparison  of  Exact  and 

Diagonal 

Buehler’s  Values, 

and  aJn, 

Respectively#  with 

the  Exact 

and  Diagonal 

Values 

UD 

a  and 

a,  ,  Respectively, 

Corresponding  to  g(x) 

-ng 

ag  “■  ■ 

^1 

^2 

anB 

adB 

a 

ng 

a  , 

5 

5 

60.7 

60.70 

60.70 

60.70 

5 

4 

51.8 

51.89 

51.89 

51.89 

5 

3 

41.2 

41.21 

41.22 

41.21 

5 

2 

31.9 

31.91 

31.91 

31.90 

5 

1 

23.3 

23.34 

23.34 

23.34 

5 

0 

12.3 

12.32 

12.32 

12.32 

4 

4 

44.3 

44.40 

44.40 

44.40 

4 

3 

35.7 

35.73 

35.74 

35.73 

4 

2 

27.2 

27.23 

27.23 

27.23 

4 

1 

18.8 

18.77 

18.77 

18.76 

4 

0 

9.05 

9.05 

9.05 

9.05 

3 

3 

28.9 

28.89 

28.89 

28.89 

3 

2 

22.0 

22.04 

22.04 

22.03 
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